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Abstract 

This paper presents asymptotic properties of the maximum pseudo-likehhood 
estimator of a vector d parameterizing a stationary Gibbs point process. Sufficient 
conditions, expressed in terms of the local energy function defining a Gibbs point 
process, to establish strong consistency and asymptotic normality results of this 

1 -^ ' estimator depending on a single realization, are presented. These results are general 

C^ , enough to no longer require the local stability and the linearity in terms of the 

parameters of the local energy function. We consider characteristic examples of such 
models, the Lennard- Jones and the finite range Lennard- Jones models. We show 
that the different assumptions ensuring the consistency are satisfied for both models 
whereas the assumptions ensuring the asymptotic normality are fulfilled only for the 

■^^ ' finite range Lennard- Jones model. 
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5 ■ 1 Introduction 
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These last years, much attention has been paid t o spatia l poin t pattern data, and especially 
k>( \ to models and methodologies fo r fittin g them, see MoUer ( 20081) for a recent overview of this 

5_j ■ to pic andlbaley and Vere-Joned (|l988l ). IStovan et all (|l987n iMaller and Waaeepetersenl (|2003l ) 



Cd \ or llllian et al.l ( 20081 ) for more general information. For spatial point pattern data, the refer- 



ence model is the Poisson point process modelling a random configuration of points with no 
interaction between points. In particular, this leads to the independence of any two random 
sub-configurations lying in two non-overlapping domains. A way to introduce dependence is 
to consider the class of Gibbs models. In a bounded domain, a Gibbs point process is defined 
through its probability measure having a Radon-Nykodym derivative with respect to a Poisson 
point process measure proportional to e^^^^' where y{}p) corresponds to the energy function 
(i.e. a cost function expressed in terms of interactions) of the configuration of points (p. The def- 
inition of Gibbs models in R'' is essential when dealing with asymptotic properties of estimators 



based on a point process observed in a domain aimed at converging towards R''. The extension 
of this definition is not so straightforward. The probabihty measure of a Gibbs point process in 
Mf^ has to be defined by specif ying its conditi onal density (indirectly expressed in terms of the 
energy function V{ip)), see e.g. IPrestonI (|l976r) or Section [2] for more details. 

The class of Gibbs point processes is extremely rich. The energy function can penalize 
points, pairs or triplets of points (see e.g. iBaddelev and Turned (|200Cl| )). More sophisticated 
models can also be obtained by considering interactions based o n the Delaunay or the fc— nearest 
neighbor graphs (Bertin ct al. ( 1999b, c)) , Voron o 'i tesse llations (jDereudre and Lavancien pOOQ)) 
or random sets (JKendall et all (|l999l) . lDereudrel (|2009l) '). 

Following the definition of a parametric Gibbs point process, the natural question of efficiently 
estimating the parameters arises. Many proposals have tried to estimate the energy function from 
an availa ble point pattern data. The most w ell-known method is the use of the likelihood function. 



see e.g. iMciUer and Waageoetersenl ( 2003 ) and the references therein. The main drawback of 



this approach is that the likelihood function contains an unknown scaling factor whose value 
depends on the parameters. This parametric normalizing constant is difficult to calculate from a 
practical point of view. From a theoretical one, it also makes asymptotic results more complicated 
to obtain. An alternative approach relies on the use of the pseud o-likelihood f uncti on. The 
idea originated from iBesaa (|l974[ ) in the study of lattice processes. iBesag et al. ( 19821) f urther 
considered this method for pairwise interaction point processes, and IJensen and Molleii (|l991l ) 
extended the definition of the pseudo-likelihood function to the general class of marked Gibbs 
point processes. The construction of the pseudo-likelihood function is based on the conditional 
densities which spare the computation of the scaling factor. 

Our paper deals with asymptotic properties of the maximum pseudo- likelihood estimator. In 
order to underline our theoretical improvements, let us discuss the two main different papers 
discussing this topic: 



• In iBilliot et al.l (J2008[l , we obtain consistency and asymptotic normality for exponential 
family models of Gibbs point processes, that is, on models with energy functions that 
are linear in terms of the parameters. Moreover, we concentrate on models such that 
the local energy function is local and stable. The locality of the local energy expresses 
that the energy to insert a point x into ip, that is, F(x|iy9) — V{(p U x) — V{ip), depends 
only on the points of (p falling into some ball with a fixed radius whereas the stability of 
the local energy (property referred as the local stability) ass erts that V( x\(fi) is bounded 
from b elow by a finite negative c onstant. The paper Billiot et al.l ( 20081) extends several 
papers (jJensen and MoUeiJ (|l99ll) . I Jensen and KiinschI (J1994I) ) and includes a large class of 
examples of practical interest: area-interaction point process, Multi-Strauss marked point 
process based on the complete graph or the fc-nearest-neighbors graph, or the Geyer's triplet 
point process to name a few. 

• Another work has been undertaken by Mase. The consistency for non necessarily stable 
local energy functions (actually for superstable and lower regular ones introduced bv lRuelld 
(jl97Gi) ) is obtained in Mase (1995i) for specific models with only two parameters -the chem- 
ical potential a nd the invers e temperature- which can be viewed as particular exponential 
family models. iMasd (J2000l) extended his work to the context of marked point processes 



and provided asymptotic normality by adding the assumption of finite range. 

Based on this literature, the main goal of this paper is to derive asymptotic properties similar 
to the ones presented before (consistency and asymptotic normality) but in a more general frame- 
work. We provide asymptotic results for general Gibbs point processes with non (necessarily) 
linear and non (necessarily) stable local energy functions. The characteristic example we have in 



mind is the Lennard-Jones model. This model, from statistical physics, is a stationary pairwise 
interaction Gibbs point process where the local energy to insert a point x into a configuration tp 
is parameterized as follows: for 6 — {9i, 62, O3) € M.^ with 6*2, 6*3 > 



V^'{x\^;e) 




Let us notice that iMasd (|1995l ) could only propose the estimation of 9i and 02 with known 63. 
The Lennard-Jones model is of great interest from several points of view. From a physical point 
of view, this model arises when theoretically modelling a pair of neutral atoms or molecules 
subject to two distinct forces in the limit of large separation and small separation: an attractive 
force at long ranges (van der Waals force, or dispersion force) and a repulsive force at short 
ranges (the result of overlapping electron orbitals, referred to as a Pauli repulsion from the 
Pauli exclusion principle). In this literature, the parameters 02 and 63 are often referred to as 
the depth potential and the (finite) distance at which the interparticle potential is zero. From 
a probabilistic point of view, this mo del con s titute s the main example of superstable, regular 



and lower regular energies studied in iRuelld ( 19701) where the author proves the existence of 



ergodic measures for such models. Fina lly, from a statistical point of view, this mode l has b een 
considered by several authors, see e.g. lOgata and Tanemural (|l981[ ). iGoulard et al.' ("19 961) for 



fitting spatial point patterns arising in forestry. In particular, let us note that, in Goular d et al. 



()1996[ ). the model is fitted by using the maximum pseudo-likelihood method. As the authors do 



not endeavour to justify the theoretical performances of the procedure, the result proposed in 
Section H] of this paper fills this gap. 

The rest of the paper is organized as follows. Section [5] introduces some background and no- 
tation on Gibbs point processes (general definitions, examples). The maximum pseudo-likelihood 
method and asymptotic results of the derived estimator are proposed in Section [3] For general 
Gibbs point processes, sufficient conditions, expressed in terms of the local energy function to 
establish strong consistency and asymptotic normality results of this estimator are presented. 
While no general condition on the model is assumed to obtain the consistency, the characteristic 
finite range of the local energy function is required to establish the asymptotic normality. For 
the sake of simplicity. Section [3] (and the resulting proofs) would concen trate on non-marke d 
Gibbs point processes. However, as we have shown in our seminal paper iBilliot et al.l ( 20081) . 



no real mathematical difficulty occurs with the introduction of marks. Section 0] focuses on the 
Lennard-Jones model. We show that the general assumptions described in Section [3] are fulfilled 
for this model. Proofs have been postponed until Section [5] 

2 Background and notation 

For the sake of simplicity, we consider Gibbs point processes in dimension d = 2. 

2.1 General notation, configuration space 

Subregions of M^ will typically be denoted by A or A and will always be assumed to be Borel with 
positive Lebesgue measure. We write A (S M^ if A is bounded. A"^ denotes the complementary 
set of A inside R^. The notation |.| will be used without ambiguity for different kind of objects. 
For a countable set J', \J'\ represents the number of elements belonging to J'; For A g E^, |A| 
is the volume of A; For a vector x S M^, \x\ corresponds to its uniform norm while ||x|| is simply 
its euclidean norm. For all x £ IR'^,/3 > and i E 1?, let B{x,p) := {y G M^, |j/ — 2^1 < p} and 
B{i,p) ■.= B{i,p)r\I?. 



A configuration is a subset <f of R^ which is locally finite in that Lpp^ :— (/? n A has finite 
cardinality N\{ip) := |(^a| for all A (s R^. The space fl of all configurations is equipped with 
the (T- algebra T that is generated by the counting variables Na{(p) with A (e M^. Finally, let 
T — {tx)x^r2 be the shift group, where t^^ : il — > J7 is the translation by the vector —x G R^. 

2.2 Gibbs point processes 

Our results will be expressed for general stationary Gibbs point processes. Since we are interested 
in asymptotic properties, we have to consider these point processes acting on the infinite volume 
M^. Let us briefly recall their definition. 

A point process $ is a il-valued random variable, with probability distribution P on (fi, J^). 
The most prominent point process is the (homogeneous) Poisson process with intensity z > 0. 
Recall that its probability measure tt^ is the unique probability measure on (fl,J-') such that 
the following holds for A d R^: (i) Na is Poisson distributed with parameter z|A|, and (m) 
conditionally to N/^ — n, the n points in A are independent with uniform distribution on A, for 
each interger n > 1. For A d R^, let us denote by n^ the marginal probability measure in A of 
the Poisson process with intensity z. 

Let 9 GM.P (for somep > 1). For any A (e R^, let us consider the parametric function Va(.; 9) 
from ri into R U {+oo}. From a physical point of view, VA^ip; 6) is the energy of ipA in A given 
the outside configuration (p\c . 

In this article, we focus on stationary point processes on R^, i.e. with T-invariant probability 
measure. For any A d R^, we therefore consider Va(.;0)) to be T-invariant, i.e. Va{tx(p;6) — 
VA(<y5; 9) for any x e R^. Furthermore, we assume that the family of energies is hereditary, which 
means that for any A isR"^, ip € ft, and x € A: Va{(p; 9)) = +oo => Va,(' P U ix}\ 9)) = + 00. 

In such a context, a Gibbs measure is usually defined as follows (see IPrestonI (|l976[) ). 



Definition 1 A probability measure Pg on fl is a Gibbs measure for the family of energies 
(Va(.; 6))a^r2 if for every A d R^, for Pa-almost every outside configuration ^pA", the law of Pa 
given ipA" admits the following density with respect to n^ : 

/a((Pa|^a.;0) = ^ .. e-^^(^^^), 

ZA((y3A<=;0) 

where Za(93a<^;^) '■— /q e~^^^''''^'~''^'^''' 'TTj^{dipA) is called the partition function. 

Without loss of generality, the intensity of the Poisson process, z is fixed to 1 and we simply write 
TT and TTA in place of tt^ and Trj^. In the previous definition, we implicitly assume the consistency 
of the family (/a(.|.; 0))AeR2: for any A C A (s R^ 

r , I ,g\_ fA{'PA^VA\A\'PA-=;0) _ fA{'PA^'PA\A\'PA-;d) 

/a(<^A\aI'^A<=;^) Jn^fAilpAUPA\A\^A'=;0)'^A{d-)jJA)' 

A sufficient condition to directly fulfill this basic ingredient is to assume the compatibility of the 
family {Va{-)) a<eR^ ■ for every A C A g R^, the function (p -^ Va{'p;9) — Va{'p;9) from fl into 
R U {+cxd} is measurable and only depends on ipA"- 

The existence of a Gibbs measure on ft whic h satisfi e s the se conditional sp e cifications is 
a diffic u lt issue. We ref e r the interested reade r to iRuelld (|l969l) : IPrestonI (|1976[ ): iBertin et al. 



(|1999al ): lDereudrel(J2005l ): lDereudre et al.l pOld ) for the technical and mathematical development 



of the existence problem. The minimal assumption of our paper is then: 



[Mod-E] : Our data consist in the realization of a point process $ with Gibbs measure Pq* , 
where 0* e 0, is a compact subset of K.^ and, for any € 0, there exists a stationary 
Gibbs measure Pg for the family {V\{.;6j)\^]g^2. 

In the rest of this paper, the reader has mainly to keep in mind the concept of local energy 
defined as the energy required to insert a point x into the configuration ip and expressed for any 
A 9 a; by 

V{x\^;e):=VAi^U{x})~VAiip). 

From the compatibility of the family of energies, the local energy does not depend on A. 
Our asymptotic normality result will require the following locality property assumption. 

[Mod-L] : There exists D >0 such that for all f e ft 

V{0\^;e) = V{0\^Bio,D)\O)- 

2.3 Example : Lennard-Jones models 

Let us present the main example studied in this paper. We call LJ-type model the stationary 
pairwise interaction point process defined for some D s]0, +oo] by 

X2£ipK<i 

and 

As a direct consequence, the local energy function is expressed as 

F^^ {x\^- 6) := 9, + i/^^ {x\^- 6) with i/^^ {x\^- 6) :- ^ g"^' {\\x - y\\-e). 

where = (01,^2, 6*3) e M x (R+)^. The cases D = +00 and D < +00 respectively correpond to 
the Lennard-Jones model (briefly presented in the introduction) and the Lennard-Jones model 
wit h finite range . 

iRuelld (J197C1I) has proved the existence of an ergodic measure for superstable, regular and 
lower regular potentials. The Lennard-Jones model (including the finite range one) is known to 
be the characteristic example of such a family of models for which Ruelle managed to prove the 
existence of ergodic measures for any S R x (M+)^. In order to ensure [Mod-E], it is required 
to assume that 6*2,^3 > 0. Finally, [Mod-L] is satisfied for the LJ-type model with D < -|-oo 
since the parameter D corresponds for pairwise interaction point processes to the range of the 
Gibbs point process. 

3 Asymptotic results of the Maximum pseudo-likelihood 
estimator 

3.1 Maximum pseudo-likelihood method 



The idea of maximum pseudo-likelihood is due to iBesasJ (|1974I ) who first introduced the concept 



for Markov random fields in order to avoid the normalizing constant. This work was then widely 



extended and I Jensen and M0lleij ( 19911 ) (Theorem 2.2) obtained a general expression for Gibbs 



point processes. Using our notation and up to a scalar factor, the pseudo-likelihood defined for 
a configuration (p and a domain of observation A is denoted by PLa (¥'; ^) and given by 

PLa {^- 9) = cxp (- f e-n-\^-'^)dx) n e'Vi.i^vfi)^ (^^ 

It is more convenient to define and work with the log-pseudo-likelihood, denoted by LPL/^ {ip; 6) 
LPLa {ip; 6) = - j e-^(^''^'^)dx - ^V {x\ip \x;0) . (2) 

The point process is assumed to be observed in a domain A„ (B D — yj^^K^^ix, D) for some 

D < +IX). For the asymptotic normality result, it is also assumed that D > D and that A„ C M.'^ 

can be decomposed into Ui^i^Ai where /„ = B (0,n) and for i G Z^, A^ = Ai{D) is the square 

centered at i with side-length D. As a consequence, as n — >■ -l-oo, A„ — >■ R^ such that |A„| — 7> -l-cx) 

, |5A„| ^ 
and , . , -> 0. 

|A„| 

Define for any configuration ip, Un{(p',0) = —j^—rLPL\^{ip;9). The maximum pseudo- 
likelihood estimate (MPLE), denoted by 6n{(p), is then defined by 

6n{(p) = arg max LPL\ (ip; 6) = arg min C/„ {p; 6) . 
Be® " 6e@ 

The following basic notation are introduced: for j', fc = 1, . . . ,p and A d R^ 

• Gradient vector of f/„: U';^\p>] 9) := -\An\-^LPL'-ll (^; 6) where 



Hessian matrix of [/„: t/l^^^'; ^) := -\A„\-'^ LPL ^l^ {p; 6) 



-A„ 
32 1 



(^i^'(^=^)),, - /,(J;]^(^I^^^)-|^(-I^=^)|^(^I^=^))^-^^"^^'^^- 



Y^ i:r\p\x-e) — {x\p^\x;9). 






Finally, note that from the decomposition of the observation domain A„, one has 

C/W(^;0) = -|A„ri Y. LPL^l] {v;0) and U^\p; 6) = -\K\-'Y.I^^Z ^^.0) 



iein iein 



3.2 Consistency of the MPLE 

The assumption [C] gathers the following four assumptions: 
[CI] For aU 9e®, 

|^g-y(ol*;0)^ < +00 and E {\V (0|$; e)\ e-^(ol*^^*)) < +oo. 



E\ 



[C2] Idcntifiability condition : there exists Ai, . . . , A^, £ > p events in il such that: 
— the £ events Ai are disjoint and satisfy PQ*{Bi) > 



— for aU (ipi, . . . , ipe) £ Ai x ■ ■ ■ x Ai 



D{0\ip,;e) = 
i = l...J 



= 0" 



where D{0\ip,; 0) := V (0|(^,; 0)-V (0|(^,; 0*) 
[C3] The function C/„((p; •) is continuous for Pn*— a.e. ip. 
[C4] For all (f ^Vl,V {Q\f] 0) is continuously diffcrentiable in and for all j = 1, . . . ,p 

2\ 



E max 
\0e& 



dV 

de., 



\^\e] 



,-y(o|*;6») 



< +00. 



Theorem 1 Under the assumptions [Mod-Ej and [C], for Pa* — almost every ip, the maximum 
pseudo-likelihood estimate 0n{^) converges towards 0* as n tends to infinity. 

3.3 Asymptotic normality of the MPLE 

For establishing the asymptotic normality of the MPLE we need to assume the four additional 
following assumptions: 

[Nl] For ah ip e n, V{O\'p;0) is diffcrentiable in = 0* . For all fc = 1, . . . , 3 and for all 
Ai, . . . , Afe, k positive integers such that J2i=i ^j — 3 and for A (e M^ 



Afc 



n 



dV 



^(o*>;r 



Ai 



-v{{-u...,o.,m0*)d^^^_a^\ <+c 



[N2] There exists a neighbourhood V{0*) of 0* such that for all ip E fl, V{O\lp;0) is twice 
continuously diffcrentiable in e V and, for all j, fc = 1, . . . ,p and e V{0*), 



E 

and 



d^v 



dOjdek 



{m-.o) 



-v{mfi] 



< +00, E 



d^V 



dOide 



jUOk 



m^;0) 



-V(O\<i>:0* 



E 



dV 



(O|$;0) 



-y(o|*;0) 



< +00. 



< +00, 



[N3] There exists Ai, . . . ,Ai, £ > p events in fl such that: 

— the £ events Ai are disjoint and satisfy Pg*{Ai) > 

— for all (tpi, . . . ,(pi) £ Ai X ■ ■ ■ X Ai the {£,p) matrix with entries ^^ {O\ipi;0*) is 
injective. 

[N4] There exists Aq, . . . , A^, £ > p disjoint sub-events oi il := {ip <E Q. : pa, = 0, 1 < |«| < 2} 
such that 



- forj = 0,...,^, P0*(Aj)>O. 

— for all (ipo, . . . , (fie) £ Aq x ■ ■ ■ x Ag the {£,p) matrix with entries I LPLj- {ipf, 6* 
i^LPL^-^ ((y9o; 0*)j . is injective, with A := UjeB(o,i)- 

The assumptions [N3] and [N4] will ensure (see Section [5] for more details) that the matrices 
U}-^\e*) and ^0,9*) respectively defined by 



([/(2)(r ))^ ^ :^ E (^gl (01$; r) gl (01$; 6*) e-^(°l*^^*)) 



(3) 



and 

s(i5, r ) = 5-2 ^ E (lpl'II ($; e*) iP£i\' ($; 6>*f ") , (4) 

ie]B(o,i) ^ '^ 

are definite positive. 

Observe that, when the ener gy function is linea r, the expressions of the assumptions [Nl] 



and [N2] arc clearly simpler (see lBilliot et al.l (|2008fl ) and that [C2] and [N3] are similar. 



Theorem 2 Under the assumptions [Mod], [C], [Nl], [N2] and [N3], we have the following 
convergence in distribution as n ^ +00 

|A„r/2[j(2)(g*) (g„($)-r)^AA(o,S(5,r)), (5) 

where ^{D,9*) is defined by Q. In addition under the assumption [N4] 

|A„|i/2s„(<i>;g„(<i>))-V2t/(^2)(^.g^^(^)) (0„(<i>)-r) ^AA(0,J^), (6) 

where for some 9 and any configuration ip, the matrix '^^{Lp; 9) is defined by 

S„(^;0) = |A„r'E E LPL^^^i^;9)LPL^l]i^;9f. (7) 

iei„ jeB(ia)n/„ 

In the following the assumption [N] will stand for the assumptions [Nl], [N2], [N3] and [N4]. 

4 Applications to the LJ-type model 

This section focuses on the LJ-type model presented in Section 12.31 and aims at proving the 
following result. 

Proposition 3 

(i) Theorem]^ holds for the LJ-type model (with D g]0, +00]^, that is for the Lennard- Jones and 

the finite-range Lennard- Jones model. 

(ii) Theorem\^ holds only for the finite-range Lennard- Jones model. 

The proof of Proposition [3] consists in verifying Assumptions [C] for the LJ-type model 
and [N] only for the finite range Lennard-Jones model. In the following, we will deal with two 
types of assumptions: 

• Integrabilility type assumptions, i.e. Assumptions [CI], [C4], [Nl] and [N2]. 



• Identifiability type assumptions, i.e. Assumptions [C2], [N3] and [N4]. 

Note that [C3] is obvious since g^^{r,) is continuous. For the integrabihty type assumptions, 
the following Lemma will be widely used. 

Lemma 4 Let ^ be a stationary pairwise interaction Gibbs point process assumed to be super- 
stable, regular and lower regular. For i — 1,2, define Hi {x\ip) — J2y£ip 9ii\\^ ~ 2/11) '"^^^ 9i ^ 
continuous function. Assume that there exists e > such that there exists a positive and decreas- 
ing function g{-) such thatg^{r) := g2{i') — £\gi(r)\ > —g{r) for allr > and J^^ rg{r)dr < +oo. 
Then for all k > 0, 



£;(|i7i(0|$)|'=e-^^("l*)j <+oo. 
Proof. For all finite configuration ip 

< c(£,fc)e-(^^("l'^)-'^^i(°l^», withc(e,fc)=f 

< c(e,fc)e-^^(ol^), 
where 



k^' 



xGip 



Now, the assumptions ensure that Qp, is l ower r egular in the Ruelle sense. We may now apply the 
same argument as in Lemma 3 of lMasd ( 1995[) to prove the integrabihty of the random variable 



Before verifying the different assumptions, let us denote by 

6*,;"^ := inf 6l„ 6*^ p := sup 6l„ 6*'"^ ;= min(6l!,"^ e'!,"^ ) and 6I""p := max(6'™P,6'^^P). 
6/e0 6>e© 

Since is a compact set of M x (]0, +oo[)2, then 6*'"^ > and 6'™p < +oo. 

4.1 Assumptions [C] 

4.1.1 Assumption [CI] 

The first part is a direct application of Lemma |31 For the second part, one has to prove that for 
all 6> e 

Let ge{r) = .g^"'(r; 0*) - e|g^"'(r; 0)|. We have 

m { .i/^ ,a ^ j if r < ^3 



m f 0*\ 0* 

which satisfies the assumptions of Lemma S] as soon as e < (#1 #, that is, as soon as 

(£)■' 



£ < 



4.1.2 Assumption [C2] 

Let us denote for n > I, Cn = B{0,n) \ B{0,n — 1) and define for m,n > 1 the foUowing 
configuration sets 

Um,n = {ip e il : \(pcj < m\Cn\} 

In order to prove [C2], we need tire following Lemma. 
Lemma 5 Let R G W^ , 9 £ & and (p G Um, let us denote by 

Z{ip,R;9):= ^ .9"(||x||; 0), 

then for all S > there exists Rq such that for all R > Rq, \Z(ip, R;6)\ < S. 
Proof. 

zi^,R;e)^\ J2 ff"(INI;^)l < E E k"(ll^ll;e)| 

< Y. I^cjx sup |g"(||x||;r)|. 

There exists a constant k = k{R) such that for all n > [i?] , sup^g(^^ |g^''(||x||; 0*)| < kn~^. 
Therefore, 

I Y. 9"mm\<km Y |C„|xn-6 = 0J ^ n 

which leads to the result since the previous series is convergent. ■ 

Let 9 E &\9* and consider the following configuration sets defined for k > 1 and for 77 small 
enough by 

Ao = {ipen:\ipr\BiO,D)\ = 0} (8) 

Akiv) = {^en:\^nBiO,D)\ = \pnBiiO,Dk-'/''),r,)\ = l}, (9) 

where D is any positive real for the Lennard- Jones model and corresponds to the range of the 
function g^'' {■) for the finite range Lennard- Jones model. There exists m > 1 such that for all 
77 > and for /c = 2, 4 

Pq* [Ao n C/„,) > and Pq* [Ak (77) n t/„) > 0. 

Now, let (^0 € ^0 n Um, <P2 e A2{vi) n Um and (^4 e Ai{ri) n [/,„. First, 

i?(0|(/7o; 0) = 01 - 0^ + Ziip^.D- 9) - Z((^o, ^; 9*) - 0. 

For the Lennard- Jones model, according to Lemma [S] one has, for D large enough, 

|Z(</7o, D- 9) - Z((^o, D- 0*)\<\ \0i ' e\\ . 
10 



Hence for 77 small enough, and for both models 

= \Dio\ipo;e)\ 

> \e,-ei\-\ziipo,D;e)-zi^o,D;e*)\ 



> 



1, 



'1 - t'lh 



which leads to 9i = O*. Moreover, 



D{0\ip2;e) = 402 2 



12 /fl N 6> 



D 



-V2 



D 



- 46*; 2 



D 
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-^'i 



-/2((^2) + ^(<P2,Ae)-^(<y52,Ae*) 



iP(0|,.;.) = 4.2 (^4 ^^j -2^^jj-4.,*(^4^^J -2^^ 

+/4((^4) + z {^4,D; e)-z ((^4, i?; e*) , 

where for any (pk G Ak{ri) (k ~ 2,4), there exists a positive function fkii]) converging towards 
zero as ry — > such that \fk{^k)\ is bounded by fkiv)- Now, we have 

4(2 - 2^/2) 



2D{0\ip2;e)-D{0\ipr,e) = 



with 



1)6 



0291 - ep;'' + 2/(^2) - /4(^4) + z'{^2, ^4, D; e, e*) 







Z'{ip2,ipi, D; e, e*) := 2 iZ{ip2,D; 6) - Z{ip2,D; 9*)) - (Z(^4, D; 6) - Z(^4, D- 9*)) 
For rj small enough, we have, for any ipk G ^fc(?/) (fc — 2,4), 



|2/((P2)-/4(<^4)|< 2/2(7/) + /4(r,)< 



4(2 - 2\/2) 



£)6 



72^3 - f2C'3 I 



For the finite range Lennard- Jones model, Z'{(p2, (p4, D; 9, 9*) = 0. For the Lennard- Jones model, 
according to Lemma [5l one has for D large enough 



\Z'iif2,^4,D;9,9*)\< 

Hence for 77 small enough, and for both models 

4(2 - 2V2) 



4(2 - 2\/2) 



D6 



72^3 - t'2t'3 I 



> 



> 



4(2 - 2^/2) 



1 



4(2 - 2\/2) 



(^2^3' - ^2^3') + 2/(^2) - /4(^4) + Z'{^2.Vi,D- 9, 9*) 

\e2el-e*ef\ - 12/(^2) - /4(^4)| - \z'{^2,'P4,d-9,9^)\ 



D6 



^2^3 - (^2^3 I 



leading to 620^ — 6*2^3®. By considering the combination V2D{0\ip2;9) — D{0\ipi;9) and using 
similar arguments as previously, one obtains: ^2^3^ = ^2^3^^- By computing the ratio of the two 
last equations, one obtains 63 = 9^ and then 62 = ^2- 
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4.1.3 Assumption [C4] 

For all ifi £ Q and for any e 0, V^'^ {0\ip; 9) is clearly differentiable in 9. First, note that [C4] 
is trivial for j — 1- For j = 2, 3, let us define: 



Xjiv;0) :- %-{o\^;0] 



^-V^''{0\^:9) 



Our aim will be to prove that for j — 2,3 and for all fc > 

J5(inaxXj(i>;6l)M < +00. (10) 

In particular, the Assumption [C4] corresponds to PH)) with k = 2. Let us notice that for all 
ip eft and for all € 



with for some r > 0, g (r) :— iO ( -^-p^ ^ ^^ I • Let us also underline that for j = 2, 3 



BnLJ I 4 ^^rW S^ if ■?' = 2, 

^^(r;6>) >gj"f(r) with g'^ffr) := <^ V. ,,,,ii ^ ^ 



Therefore, by defining Vj'^^{0\ip) := J2xeip9T^i\\^\\)y the result ((TO)) will be ensured by proving 



According to Lemma HJ in order to prove this, let us denote by gj,e{) the function defined for 
j = 2, 3, for some e > and for r > by 5j,e(r) — ]jT^{r) — e |g'"^(?')| • On the one hand, one has 

* I rl2 r6 I 11 ' ^ gaup , 

which satisfies the assumptions of Lemma H] as soon as e < 6™^ . On the other hand 



4^- , v^-)"-l^K^-)" _ io-n^-^io-n^] if , < ^^f-""^ ^'^ 



'"^^'^ ^ 4,.. / (^-rwKir _ (^-)-.e.(r^.)A .^ ^ ^ /^ (^^ ^ ^/^ 



which satisfies the assumptions of Lemma H] as soon as e < 6''"^/12, which ends the proof. 
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4.2 Assumptions [N] 

4.2.1 Assumption [Nl] 

Let us present two auxiliary lemmas. 

Lemma 6 Let (p be the realization of a stationary pairwise interaction point process with local 
energy function defined by 

vix\ip;e)^ei + H{x\ip;e) with H{xW;e)^Y.9{\\v-A\:0). 

Let K < +00 and let xi, . . . , xk G M."^ \ip, Xi y^ Xj for i,j — l,...,K (where K < +ooj, then 

K 
H{{xi,...,XK}W-e) = ^i/(x-fc|<^;0) + i/({xi,...,x;f};0) 

fc=i 

K 

V{{xi,...,XK}W-e) = ^V^(xfc|^;0)+i/({a:i,...,a:K};0) 

fc=i 

This result comes from the definition of the local energy. 

Lemma 7 Using the same notation and under the same assumptions of Lemma\^ assume that 
there exists gmin such that for all r > and any 6 (z &, g{r; 6) > gmin, then 

fc=l 

Proof. The proof is immediate since 

H {{xi,...,xk};0) ^^g{\\xz -Xj\\;e)> g™„. 

■ 

Let k — 1, . . . , 3 and let Ai, . . . , A/c, k positive integers such that J2i=i ^i = ^ ^'^d define the 
random variable 
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Cfc 



< Cfc 



A*-' 



A*-' 



= 1 



From Lemma [3 we have 

k 

^^ 

fc 

cfciAi'=n^ 



dV_ 
dV_ 



(xj$;r) 



(xj$;r) 



-y(xi|$;0*) 



da;, 



'yi-^\'^:n]dx^...dxk 



i'-l*^"-) 



l/fc 



^y(.d*;r) ^^^^ 



■ dxk 



1^(^.1*;^ 



9y 



fc \ i/fc 



fc \ i/fc 



by using Holder's inequality and the stationarity of the process. The result is then a simple 
consequence of pUj) and Lemma 2] 

4.2.2 Assumption [N2] 

For all If ^ n, it is clear that for all 9 £ &, V{0\(p;O) is twice continuously differentiable in 9. 
According to Lemma 3] and the fact that [Nl] is satisfied, it is sufficient to prove that for all 
J, fc=l,2,3 

g2yL.J 



iO\<P;9] 



■'{0\<!>;9) 



dOjdek 

q2 LJ 

This is obvious when either j or k equals 1 and when j = k = 2 (since jog-yxii", •) — 0). Now, for 



the other cases, define for e gj^k,s{i~) '■= 9 {f", 9) — e Qg^gg {r; 9) 
92,z,e{r) = g3,2,s{r) = 



We have 



4 C^2e;^"-12693^' ot~6e9n n r < 2^/^ 



92e" + 12e9" _ e|+6ee£ 



otherwise 



„t>2 



which satisfies the assumptions of Lemma Has soon as £ < ^j|^, that is, as soon as e < ^ ^^ 
Finally, 



'82 9,^^-13260.^° eaea-SOee ;^^ ;f „ ^ /132\l/6, 



53,3,eW 



9291^ + 132601" _ e2e°+3ae9^ 



ifr<(^) 
otherwise 



which satisfies the assumptions of Lemma S] as soon as e < -j^, that is, as soon as e < 



„t^3 



(^ 

132 



4.2.3 Assumption [N3] 



Let y = (2/i,y2,2/3) G K^ and g{y,ip) := y^F^j (0|^; 6»*). Let ipo € Aq and (^^(ry) € Akiv) 
(fc — 2,4) where Aq and Ai;(r/) are defined by ([SJ and ©• Assume g{y,(pk) — ior k ~ 0,2,4. 
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Since, g{y, (po) — yi, we have yi — 0. Now 

9{y,'P4) 



4.,|2||)"-V^(| 



1261 



66i; 



4y3^2M2^T^-V2-^ +/2(y,¥'2) 



D 



D^ 



42/2 



D 



12 



4(13 _2M^ 



D 



4y3e2 4 



126i; 
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£> 



12 



,6^ 



+ f4{y,V4:), 



where for any ipk G ^k{'>]) (^ = 2,4), there exists a positive function fkiVi'n) converging towards 
zero as ?7 — > such that \,fk{y, <Pfc)| is bounded by ,fk{y, v)- Now, we have 

My, V2) - g{y, Pi) = 4(2 - 2V2))^ i9;y2 + 6e*y3) + 2h{v, ^2) - h{y, Pi) = 0. 
For 77 sniaU enough, we have, for any ipk G Akiv) (^ = 2,4), 

1 



|2/(y, ^2) - /4(y, ^4)1 < 21/2(2/, r;)| + |/4(y,r,)| < 



4(2 -2V2)-^(03*?y2 + 60^2/3) 



Hence for 77 small enough. 



0= \2g{y,ip2)-g{v,Pi)\ > - 



4(2 -2V2)-^(0*y2 + 60^2/3) 



leading to the equation 9t,y2 + 6^22/3 = 0- By considering the linear combination \/2g{y,p2) — 
g{y, Pi), we may obtain the equation 0^y2 + 12^22/3 — with similar arguments. Both equations 
lead to 2/2=2/3 = 0- 

4.2.4 Assumption [N4] 

The assumption [N4] may be rewritten for all fc = 1, • • • ,£ and for all ipk G Ak and ipQ (L Aq: 



Vy e m3, y^ [I^PI^-K ^^^^ ^*) " ^^^4 (^o^ ^*)) = ?/^(-^(^/c; ^ ) - R{Pk;e*)) = Oj ^ y = 0. 
where for any configuration (^ G il and (^0 G Aq 

L{^]9*) := / Fi'J(a;|(^;r)e-^^'("l^'^*)dx- / V^'j (a;|</Jo; 0*) e"^'''("l^'''^*)(ia; 

Ja ' J K 

Ri^;0*) := Y._'^Lli^\f\^;()*)- 5]_Fi'](x|^o\x;r). 

Concerning this assumption, we choose (po £ Aq ~ {p £ ft : (^Aq =0}- Let y G K^ then 

f_y'V^l]{x\ipo;e*)e-'''-'i-\'^»-'(^*)dx^y,e-'>*^\A\ and ^ y^V « (x|^o \ x; r) = 0. 

Ja t- 

icGi^onA 

Consider the following configuration set, defined for ry, e > 0, by 

^2(17, e) = {<p G il : </?Ao = {zi,Z2} where zi G S(0,?7),Z2 G B{{0,2ri + e),r])} . 
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Note that for zi G B(0,7?),Z2 € 6((0,2r/ + e), 77), e < \\z2 - zi\\ < e + 477. Let ip2 e Miv^^) and 
a; S A, then one may prove that for j = 2, 3 

y"(x|^2;r) = e* + 25"(||x||;r) + /(x,77,e) 

-^^{xW^-X) - 24^(||x||;r) + /,(:r,r;,£) 



90, 



where f{x,ri,e) and fj{x,'r],e) are such that 



hm f(x,ri,e)= hm f,(a;,77,e) = 0. 

On the one hand, one may prove that there exists a function /^(y, ?], e) such that 
hm(^„)^(o,o) /L(y,?7,e) = and such that 

y^i(^2; r ) ^1/1- yie-^'i |A| + /^(y, ij, e) 
where 

7:=j^MI|a;||;e*)e-'^-'^"^"""^^*^rf^ and h{r-e*):= A, 2^(r; 0*), 2^(r; r) 
On the other hand, there exists a function fR{y, 7], e) such that hm^_i.o fniy^ ?y, e) =0 



/il((^2; 0*)=2yi+ 22/24 ' ''''' '^^ 



^y^^o*2{j^-'-^)+My,v,e). 



Since 



~E^ {2V2^Qf + 2y34026^35) + 2y240f ' + 2^340^120^^ 



For 77 and e chosen small enough, one may prove that 



^=y^f{L{^2;e*)-R{^2:e*))\ > 77 2y2A9f^ + 2y3W;i20;^^ 



leading to 

This means that 



2y240f ^ + 22/340J120*" = ^ 0*2/2 + I29^y3 = 0. 



(11) 



y^i2(¥.2; e*) = 2j/i - ^ (2^2403*' + 2^340^60^') + fRiy,v,£)- 
With similar arguments, we obtain that 



2?;240f + 2y3Ae*60f = ^ O^y^ + W^y:, ^ 0. 



(12) 



Equations ([TT]) and ([T^ lead to 1/2=2/3 = 0. Now consider the following configuration set defined 
for some fc > 1 and 77 > 

Mv) = {v' e ^ : '/'Ao = |<P n 6(0, 77)1 = fc} 
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and let ipk G Afe(ry). Then, one may prove that there exists a function fLiyiil) such that 
Uni^^o fhiVi'n) — ^-nd such that 

y" {L{^k;e*) - R{^k-X)) = yi j_e-'^' (e-'^-^"'(ll^ll^^*) - l) dx - ky, + /^(y, ry) = 0. 

Let us denote by Ai :— B{0, min(6'3, D)) and A2 := B{0, D) \ Ai Now let us consider two cases. 
Case 1: ^3 < D. First note that for all x G A, (/^■■'(||x||; 6*) > 0. Then, for k large enough and 
for 77 small enough, we have 



-n(e-^9^''(\\-\V0')_i\dx 



|Ai| e* 1 
< V^e^^i < - and 
k 4 



^fL{y,v) 



< 



\yi\ 



Hence for k large enough and for rj small enough, we may obtain 

= j\y'{L{^k;e*)^R{^k;e*))\ 



k 
> \yi\ 



yi 



IJ^ e-«^*(e-^-^"(ll-ll^^*)-l)dx + i/.(y,r;) 



^ I I \yi\ \yi\ I2/1I 

which leads to j/i = 0. 

Case 2: 9^ > D. First note that for ah x € A2, 



^LJf 



;©*)< 



^"(Ar) = 40*((^)"-(|) )<o. 



On the one hand, for k large enough and for 77 small enough, we may have 



lyi J^ e-^'^ (e-^9^-'i\\^wfi') _ 1) dx + IfUy, V) 



2/1 



< 



yil < ^\yi\ 



On the other hand, we have for k large enough 



y^ I e-'^J(e-^9"(ll-ll^^*)^l)da; 



\yi 
k 



> 

- k 

> 2|yi|. 
Therefore for k large enough and for rj small enough, we have 



7/1 I ^ ^ pfelSmI _ 1 

^e-^^>2| {e-^^- - 1) = \yi\^-'^-^-^ 



= i|y^(i(¥>fe;r)-il(^,;r))|>2| 



yi|- -;^\yi\ 



2 



which leads to j/i = 0. 



5 Annex: proofs of Theorems [T] and [2] 

Let us start by presenting a particular case of the Campbell Theorem combined with the Glotz 
Theorem that is widely used in our future proofs. 
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Corollary 8 Assume that the point process <& with probability measure P is stationary. Let 
A (E M'^, ifi £ fl and let g be a function satisfying g{x,ip) = g{0,Tx(p) for all x € M.'^ . Define 
f{ip) = 5(0, V5)e-^(0|'^) and assume that f G ^^(P). Then, 



E i Y. ff(^' * \ ^) ) ^ ^ (X ^^''' *)e"''^''*^^=^) = |A| ^( .9 (0, <&) e 



-\/(o|$) 



(13) 



Proof, sec Corollary 3 of iBilliot et alj (|2008D ■ 

Let us now present a version of an ergodic theorem obtained by iNguven and ZessirJ ( 19791 ) 
and widely used in this paper. Let Ag be a fixed bounded domain 



Theorem 9 (Ngu yen and ZessinI ( 1979|)) Let {Hq, G € Bt} be a family of random variables, 
which is covariant, that for all x G R^, 

Ht,,g{txV) = Hciifi), for a.e. (p 

and additive, that is for every disjoint Gi, 6*2 G Bb, 

HG1UG2 = Hgi + Hg27 a.s. 

Let X he the sub-a -algebra of T consisting of translation invariant (with probability 1) sets. 
Assume there exists a nonnegative and integrable random variable Y such that \Hg\ <Y a.s. for 
every convex G C Aq. Then, 



lim 



1 



for each regular sequence Gn — > R^- 
5.1 Proof of Theorem [1] 



i^G^ 



1 



|Ao 



■S(ffA„|I), a.s. 



Due t o the decomposition of stationary measures as a mixture of ergodic measures (see iPreston 
( 19761 )). one only needs to prove Theoreni[T]by assuming that Pa* is ergodic. From now on, Pg* 
is assumed to be ergodic. The tool used to obtain the al most su r e conv ergence is a convergence 
theorem for minimum contrast estimators established bv iGuvonl ( 1992 ). 
We proceed in three stages. 
Step 1. Convergence ofUn{^;0). 

Decompose C/„((p;0) = ^ {Hi^A^{(p) + iJ2,A„(<p)) with 

Hi,AAf)= f e-^(^\'^--^)dx and H2,aAv) = Y, V {x\ip\x;e) . 



^S*A„ 



Under the assumption [CI], one can apply Theorem!^ ( Nguven and ZessirJ ( 19791 )) to the process 
_ffi_A„- And from Corollary [51 we obtain P^*— almost surely as n — >■ +00 



|A, 



-i/i,A„($)^£; e 



'(O|*;0) 



(14) 



Now, let G C Aq, we clearly have 



x^ipG 



xeipi^, 
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Under the assumption [Mod] and from Corollary [51 we have 



J2 \V ixl'P \x;e)\\= \Ao\E {\V (0|$; 6) |e-^(°l*^^*)) < +00 



This means that for all G C Aq, there exists a random variable Y G L^jP^*) such that 
H2.g (^)\ < Y. Thus, under the assumption [CI] and from Theorem [9l ( Nguven and ZessinI 




( 19791 )) and from Corollary [51 we have P^*— almost surely 



1 „ ,,. 1 



H2,aA'^)^j^e( Y, ^(a;|$\x;0))=i;(T/(O|$;0)e-^(oi*^^*)). (15) 



We have the result by combining ([Hj) and ([T5|) : P^*— almost surely 

[/„($; 9) ^ U{e) = £;(e-^(ol*'^) + V (0|$; 6) e-^(ol*;^*)^ (16) 

Step 2. Un{-',6) is a contrast function 

Recall that £/„(•; 9) is a contrast function if there exists a function K{-, 0*) (i.e. nonnegative 
function equal to zero if and only if = 9*) such that Pg* —almost surely t/„ ($; 0) — C/„ (<&; 0"^) — > 
K{e,e*). From Step 1, we have 

K{e, r)=^(e-^(°l*'^*) (^e^(o|*;0)~v(o|$;r ) _ j^i + ^ (01$; 0) ^ y (0|<i>; 0-^) ))) . (17) 

Since the function i 1— )■ e* — (1 + t) is nonnegative and is equal to zero if and only if i — 0, 
K[e,e*) >0 and 

K{e,e*)=Q ^ e^(°l^'^)-^(°l'^'^*)- (l + y(O|(^;0)-y(O|(/j;r)) =0 
4^ i:>(0|</?;6/) := V"(0|(y9;6/)- V"(0|(y9;6/*) =0 

for Pq* — a.e. (/9. Let us consider the i events Aj (j = 1, . . . , £) defined in Assumption [C2]. The 
previous equation is at least true for ipj e Aj, which leads under Assumption [C2] to 6 = 0* . 
Therefore, K{9,9*) = Q ^ 6 = 6* . The converse is trivial. 

Before ending this step, note that the assumption [C3] asserts that for any 93, Un{if; •) and 
K{-,9*) are continuous functions. 
Step 3. Modulus of continuity. 

The modulus of continuity of the contrast process defined for all (^ G fJ and all 77 > by 

W^„((^, ri) = sup {\Uniip; 9) - Univ, 0')\ .0,0' G@,\\e 0'\\ < r;} 
is such that there exists a sequence {sk)k>i, with e^ — )• as fc — > +00 such that for all A; > 1 

P (^limsup (Wn U, l^ ^ £fc) ) = 0- (18) 

Let us start to write Wn (<y3, ■^) < Wi^n {f, ^) + W^2,n ("yS, ■^) with 

Wi,nU,l) := snpiw[j,Jip-0,9'):9,9' e&,\\0-0'\\<^\ 

:- sup fwij,J^;0,9') : 9,9' e @,\\0 - 0'\\ < l\ . 
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W2,n [ if, -^ 



and 



wi.{^-e,6') 



|A, 



n\ JK 



'V{x\v;6) _ -V{x\v:0') 



dx 



wi^^j^;e,e') := -^ ^ \v{x\^\x-e)-v{x\^\x;e') 



Let fc > 1 and let 6,6' E & such that ||0 — 0'|| < ^, then under the assumption [CI] and from 
Theorem [9] and Corollary |8l we have P^*— almost surely as n — > +oo 



w^i,A„('J';^'^') 



E 



,-V {0\<S>;6) _ -v{o\<s>;6') 



W{^J<^;e,e') — ^ £;(|l^(0|$;6»)-F(0|$;6>')|e~^(*'l*^^*) 



Under Assumption [C4], one may apply the mean value theorem in MP as follows: there exist 
^''^\ ..., l'^) € 11^=1 [min(6'j, ^j), max(0j, 61^.)] such that for all (^ € fi 

This leads, under Assumption [C4], to the following inequality 



-y(o|*;0) _ -V(O|*;0') 



r ^< 



-V{p\<t>;0) _ -y(o|*;0') 



< ^ 



(ll-^'lPE||^(0|^;.-)e-(-^"')f.) 



^ 'i' ^- 



with 71 ■.= E(Y.Um^^Q^Q\§-m:e)e-''^'>\''-(^) 
prove that 

E (\y (0|$; e)~V (0|$; 0') I e-^(oi*'^*))' < ("i 



< +00. In such a way, one may also 



ll 



with 72 := £;fe^imax0^0 ||1 (0|ci>; 0) e-^C^I^^^*) 



Hence, for all A: > 1 and for all 



6,6' E & such that ||0 — 0'|| < -^ there exists no(fc) > 1 such that for all n > no(fc), we have 

H^^^ (^; 6, 6') < ^71 and Wl^^^ (</.; 0, 0') < ^72, for Pg. - a.e. if. 

Since 71 and 72 are independent of 6 and 6 , we have for all n > no(fc) 

W„ L, i j < VKi^„ (^. ^j + W^2,« U,lj <-^ (71 + 72) := p for Pq* - a.e. <f. 



Finally, since 

limsup<^ Wni(p,T] >T 



n u{»'-(^.i)si}c u {»"(-^)^i} 



mGN n>m 



n>no(k) 



20 



for Pg*— a.e. ^p, the expected result (IT5|) is proved. 

Concl usion step. The Steps 1, 2 and 3 ensure the fact that we can apply Property 3.6 of lGuvonI 



(jl992n which asserts the almost sure convergence for minimum contrast estimators. 



5.2 Proof of Theorem M 

Step 1. Asymptotic normality of Ul-^' {^;6*) 

The aim is to prove the following convergence in distribution as n — > +cxd 

|A„|i/2 t/W($;r) ^ AA(0,S(5,r)) (19) 

where the matrix S(D, 6*) is defined by (U). 

The idea is to apply to Ull'{^;0*) a central limit theorem obtained by I Jensen and Kiinsch 



(|1994l ). Theorem 2.1. The following conditions have to be fulfilled to apply this result. For all 






(i) For aU ieZ^,E[ [LPV^[ ($; 6*)] |$aj - 



{a) For &\\i(^I?, E 



3N 



LPU^'i^-^e") 



< +00. 



(Hi) The matrix ¥ar f |A„|i/2C/^i'($; 0*)\ converges to the matrix ^{D, 9*). 
Condition (i) : From the stationarity of the process, it is sufficient to prove that 

(zPzW(cI>;r)) |<i>A5)=0. 



E 

Recall that for any configuration tp 
( 



LPL^ll{^;e*))^ = -J^ |^(x|^;r)e-^(-l'^^^*)dx + ^ |^ (x|^ \ x; ^ ^(d^)- (20) 

Denote respectively by Gi('/?) and 6*2(1^9) the first and the second right-hand term of (|20p and 
hy Ei = E (G'i($)|$A= = 'PA^)- Let us define for any ip, the measure fi^p := X^xg^j ^^- F™!^ the 
definition of Gibbs point processes, 

E2 = ^ r / 7rA„(d^A„) / M^.„(rfa;)lAo(a^)|^ (a:|^ \ x; 0*) e-^-o(^^^*). 

Since tt is a Poisson process, 

T^Ao{dVAo)f{<P) ^ T^Aoid'PAo) j TTA-{dip'^a)f{(p) 

and therefore, by introducing ijj := ip^^ U ip'^c 

E2 = ^^ r / 7r(dV) / A^^(dx)lAo(x)§^ {x\^A,yJVAi\x;e'')e-''^"^^^°'~'''-'^-^^*). 
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Now, from Campbell Theorem (applied to the Poisson measure n) 



^Ao(9'As) Jao Jn 



,e* 



de 



where tt^ stands for the reduced Palm distribution of the Poisson point process. Since from 



Slivnyak-Mecke Theorem (see e.g. lM0ller and Waagepetersenl (|2003l )). tt = tt^,, one can obtain 



-E2 



dV 



T{d%p) I dx — {x\ipAa U VA- ; 0*) e 



-Vaq (icUl/JAoUl^Ag;^*) 

^Ao(¥'A|i) Jn 

= y-^J -Ae(d^A„)/ dx|^(x|^;r)e-^(^l-^*)e-^-o(.;^*) 

= -El 
Condition {ii) : For any bounded domain A one may write for j = 1, . . . ,p 



LPL^l\<^;e*\ 



< 4 



J^^ix\<,;e*)e-H^^-^ndx 






a;6i/3A 



The assumption [Nl] ensures the integrability of the first right-hand term. For the second one, 
note that 



T, 



J2 ^ix\^\x;9*) 



< 



a;e*A ■' 

E 

Xi^Xi,X2^X3,X2^X3 
Xi.X2&ipi\.Xi^X2 



dV 



{xi\(p\xi;6*) 



dV 



ix2\^\x2;0*] 



dV 



{X3\ip\x3;e*] 



dV 



{xi\(p\xi;e*) 



dV 



{x2\>f\x2;e*) 



E 



dV 

— {x2\ip\xi]e*) 



The result is obtained by using the assumption [Nl] and iterated versions of Corollary |8l 

Condition (iii): let us start by noting that from the assumption [Mod-L], the vector LPL-^[ {ip; 6*) 

depends only on ipA^ for j e IB {i, 1). Let E,^j := E (lPL^I] ($; 6*) LPL''^^ ($; e*)\. Based 
on our definitions, we have 



r(|A„|i/2C/W(<i>;r)) - |A„riVar ( ^ iPi^ ($; r) ) 

= |A„|-i Y. ^^•^ 



ij(^l-n 



= lAnp' E E ^^.- + E E^^^ 

iei„ \je/„nB(i4) je/„nB(i,i)'^ 
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Let j G /„ n B {i, ly , since LPL ^[ ((^; 9*) is a measurable function of (/Ja^, we have by using 
condition (i): 

E [LPL^l] ($; 6") LPL^l] ($; 6»*f j = E (e [lPL^I] ($; 6/*) iPi^^] ($; e**)^ |$Af 



■(1) 



E E [LPV^[ ($; 6/*) |$a^ ] i-P^k (*; 0* 



-(1) 



Denote by /„ the fohowing set 



We now obtain 



/„ = /„n(U,ea/„IB(i,l)). 



Var(|A„|V2c/(i)(ci,;r)j = |A„ri ^ ^ i?,, 

ie/n je/„nB(i,i) 



= lAnrM E E ^^,. + E E ^^-^ 

\iei-a\ir, je/„nB(i,i) ^e/^ ie/„nB(i,i) 

Using the stationarity and the definition of the domain A„, one obtains 

l^"l"' E E ^-^■ = |A„r'|/„\/„| E Eo.j^UD.e*) asn^+oo 



and 

lAnP' E E ^^^^ 

Hence as n — >■ +oo 



Var 



(|A„r/2j7«($;r) 



< |A„| V«l E I^OjHO asn^+oo. 
jeB(o,i) 



=- iA„r^E E ^^^- 

ie/„ ie/„nB(i,i) 

"^°" i/„iiA„ri E Eo,k=uD,e*) 



(21) 



Step 2. Domination of Uj^' {^;0) in a neighborhood of 6* and convergence of [/j M^; d*) Let 
j,k ~ l,...,p, recall that I U)^'{ip\ 6) 1 is defined in a neighborhood V{9*) of 6* for any 
configuration ip by 



{u^Hv;0)) 



1 f d^V 



j-.k |A„| Jf^^ dOjOOk 



(x|(p; 6) exp (— F (x|(p; 6)) dx 



+ P^ / S- (^'^^ ^) 1^ (^1^' ^) '^^P ^"^ ^^1^' ^)) '^^ 



1 ,^ d^V , 



|A„| ^=^ 96i,a6'fe 



(22) 
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Under the assumption [Nl] and [N2], from Theorem 1^ ( Nguyen and ZessinI ( 19791 )) and from 
CoroUary[51 there exists no G N such that for all n > no 



{ui'\^;0) 



j,k 



< 2E 

+2x E 



d^V 



de.dOk 






im-.o) 



dV dV 



desOk 



\^-e) 



86^ 

^-v(a\ii-.9* 



dOk 



\^\e\ 



^-V(O|*;0) 



Note that from Theorem [51 ( Nguyen and ZessinI ( 1979| )). U)^'{-\ 0*) conyerges almost surely as 
n — ?► +00 towards U}'^'{6*) defined by ^. Note that U}'^'{6*) is a symmetric positive matrix 
since for all y eMP 

i/u^'Heny = E ('(y^y«(o|$;r))%-^(oi*^^*)') > o, 

where for j = 1, . . . ,p, ip € fl and for 9 e V{9*) (v^^^\x\ip;9*)) := §f{x\ip;9) and it is a 



definite matrix under the assumption [N3]. 

Conclusion Step Under the assumptions [Mod] and [Ident], and using Steps 1 and 2, one can 
apply a classical r esult concern ing asymptotic normality for minimum contrast estimators e.g. 
Proposition 3.7 of lGuyonI (J1992I) in order to obtain ([S])- 



It remains to prove (|6]). This may de done in two different steps. The first one consists in 
verifyi ng the positive defi niteness of the matrix S(_D,0*). The proof is strictly similar to the 
one of lBilliot et al.l ( 2008) (p. 261) except that the assumption [SDP] is now simply replaced by 
the more general one assumption [N4]. Now, the convergence in probability of S„(i';0„($)) 
towards S(I?, 9*) is obtained by applying Proposition 9 of ICoeuriollv and Lavancieii (|2010l ). 



Acknowledgements 

We are grateful to the referee and the associate editor for their comments which helped us in 
improving a previous version. 

References 

A. Baddeley and R. Turner. Practical maximum pseudolikelihood for spatial point patterns (with 
discussion). Australian and New Zealand Journal of Statistics, 42:283-322, 2000. 

E. Bcrtin, J.-M. BiUiot, and R. Drouilhet. Existence of "Nearest- Neighbour" Gibbs Point Models. 
Ann. Appl. Prohah., 31:895-909, 1999a. 

E. Bertin, J.-M. BiUiot, and R. Drouilhet. Spatial Delaunay Gibbs Point Processes. Stochastic 
Models, 15(2):181-199, 1999b. 

E. Bertin, J.-M. Billiot, and R. Drouilhet. fc-Nearest-Neighbour Gibbs Point Processes. Markov 
Processes and Related Fields, 5(2):219-234, 1999c. 

J. Besag. Spatial interaction and the statistical analysis of lattice system. J. R. Statist. Soc. Ser. 
B, 26:192-236, 1974. 



24 



J. Besag, R. Milne, and S. Zachary. Point process limits of lattice processes. Ann. Appl. Prob., 
19:210-216, 1982. 

J.-M. Billiot, J.-F. Coeurjolly, and R. Drouilhet. Maximum pseudolikelihood estimator for ex- 
ponential family models of marked Gibbs point processes. Electronic Journal of Statistics, 2: 
234-264, 2008. 

J.-F. Coeurjolly and F. Lavancier. Residuals for stationary marked Gibbs point processes, sub- 
mitted for publication, 2010. http : //hal . archives-ouvertes . f r/hal-00453102/f r/. 

D. Daley and D. Vere- Jones. An introduction to the Theory of Point Processes. Springer Verlag, 
New York, 1988. 

D. Dereudre. Gibbs Delaunay tessellations with geometric hardcore condition. J. Stat. Phys., 
121(3-4):511-515, 2005. 

D. Dereudre. The existence of quermass-interaction processes for nonlocally stable interaction 
and nonbounded convex grains. Adv. in Appl. Probab, 41(3):664-681, 2009. 

D. Dereudre and F. Lavancier. Campbell equilibrium equation and pseudo-likelihood estimation 
for non-hereditary Gibbs point processes. Bernoulli, 15(4):1368-1396, 2009. 

D. Dereudre, R. Drouilhet, and H.O. Georgii. Existence of Gibbsian point processes with 
geometry-dependent interactions, submitted, 2010. http://arxiv.org/abs/1003.2875. 

M. Goulard, A. Sarkka, and P. Grabarnik. Parameter estimation for marked Gibbs point pro- 
cesses through the maximum pseudo-likelihood method. Scandinavian Journal of Statistics, 
23(3):365-379, 1996. 

X. Guyon. Champs aleatoires sur un reseau. Masson, Paris, 1992. 

J. Illian, A. Penttinen H., and Stoyan. Statistical analysis and modelling of spatial point patterns. 
Wiley-Interscience, 2008. 

J.L. Jensen and H.R. Kiinsch. On asymptotic normality of pseudo likelihood estimates of pairwise 
interaction processes. Ann. Inst. Statist. Math., 46:475-486, 1994. 

J.L. Jensen and J. M0ller. Pseudolikelihood for exponential family models of spatial point pro- 
cesses. Ann. Appl. Probab., 1:445-461, 1991. 

W. S. Kendall, M. N. M. Van Lieshout, and A. J. Baddeley. Quermass-interaction processes 
conditions for stability. Advances on applied probability, 31:315-342, 1999. 

S. Mase. Consistency of maximum pseudo-likelihood estimator of continuous state space gibbsian 
process. Ann. Appl. Probab., 5:603-612, 1995. 

S. Mase. Marked gibbs processes and asymptotic normality of maximum pseudo-likelihhod esti- 
mators. Math. Nachr., 209:151-169, 2000. 

J. M0ller. Parametric methods for spatial point processes. Technical Report Research Report 
R-2008-04, Department of Mathematical Sciences, Aalborg University, 2008. 

J. M0ller and R. Waagepetersen. Statistical Inference and Simulation for Spatial Point Processes. 
Chapman and Hall/CRC, Boca Raton, 2003. 



25 



X.X. Nguyen and H. Zessin. Ergodic theorems for Spatial Process. Z. Wahrscheinlichkeitstheorie 
verw. Gebiete, 48:133-158, 1979. 

Y. Ogata and M. Tancmura. Estimation of interaction potentials of spatial point patterns through 
the maximum likelihood procedure. Annals of the Institute of Statistical Mathematics, 33(1): 
315-338, 1981. 

C.J. Preston. Random fields. Springer Verlag, 1976. 

D. Ruelle. Statistical Mechanics. Benjamin, New York-Amsterdam, 1969. 

D. Ruelle. Superstable interactions in classical statistical mechanics. Commun. Math. Phys., 18: 
127-159, 1970. 

D. Stoyan, W.S. Kendall, J. Mecke, and L. Ruschendorf. Stochastic geometry and its applications. 
John Wiley and Sons, Chichester, 1987. 



26 



